function[u0]=pdeic(x)

u0=[2*sin(abs(pi*x)) ; 3 * sin(abs(pi*x))]

%u0=[1;0];
%u0=[0;1];
%u0=[sin(pi*x);sin(pi*x)];

